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ABSTRACT 

This  study  contains  the  analysis  of  the  costs  associ- 
ated with  the  printing  operation  of  the  messages  received 
from  the  Naval  Communications  Fleet  Multi— Channel  Broad- 
cast System  by  a  ship  equipped  with  the  Multiplex  General 
Address  Reading  Device  (CARD) .   A  computer  simulation  is 
used  to  compare  the  average  message  delays  for  three  dif- 
ferent models  of  printers  in  varying  quantities.   Messages 
are  accounted  for  by  all  four  precedences.   Costs  of  mes- 
sage delay  and  initial  purchase  price  are  used  to  compare 
the  alternatives.   The  results  of  the  study  show  that  two 
high  speed  printers  such  as  the  TT— 624  is  the  optimal  con- 
figuration . 
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I .   INTRODUCTION  AND  STATEMENT  OF  THE  PROBLEM 

This  study  concerns  the  costs  associated  with  the 
printing  operation  of  the  messages  received  from  the 
Naval  Communications  P'leet  Multi— Channel  Broadcast  System 
by  a  ship  equipped  with  the  Multiplex  General  Address 
Reading  Device  (GARD) .   The  Multiplex  GARD  unit  performs 
the  function  of  gathering  all  messages  that  are  received 
by  the  various  channels  and  then  feeding  the  messages  to 
a  printer  or  printers.   The  GARD  unit  also  performs  the 
function  of  "weeding  out"  messages  that  are  not  addressed 
to  that  particular  ship  by  reading  the  coded  heading  and 
then  only  sending  those  messages  to  the  printer  that  are 
addressed  to  that  ship.   This  aspect  of  the  GARD  was  its 
primary  reason  for  development.   However,  the  GARD  also 
eliminates  the  need  for  exactly  as  many  printers  as  chan- 
nels that  the  ship  is  monitoring.   In  the  case  of  a  des- 
troyer, this  is  four  channels  and  four  printers.   The 
problem,  then,  becomes  a  queueing  problem  with  the  GARD 
acting  as  the  waiting  line  area  and  the  printer  or  print- 
ers acting  as  the  service  channels  in  the  model.   The  ob- 
jective is  to  minimize  cost,  where  cost  includes  a  charge 
for  message  hours  lost,  over  a  range  of  the  number  of 
printers  and  also  three  types  of  printers  that  are  cur- 
rently in  stock  or  are  being  evaluated  for  future  purchase 


The  three  types  of  printers  that  were  investigated  in 
this  study  were  (1)  the  Model  28— the  present  teletype 
used  in  most  ships  with  a  speed  of  100  words  per  minute 
(wpm) ,  (2)  the  Diablo  Model  1200  Hytype  I  Serial  printer 
with  a  speed  of  3  60  wpm,  and  (3)  the  TT— 624  produced  by 
Data  Products  with  a  speed  of  3000  wpm.   The  first  was 
looked  at  for  obvious  reasons.   The  Diablo  printer  has 
been  tested  to  some  degree  by  the  Naval  Electronics 
Laboratory  Center  (NELC) ,  in  San  Diego,  California,  and 
is  considered  a  candidate  for  future  purchase  by  the 
Navy.   The  TT— 624  has  undergone  quite  extensive  testing 
as  a  high  speed  printer  in  the  laboratory  as  well  as  at 
sea. 

The  problem  is  to  measure  in  dollars  the  various 
costs  associated  with  the  different  printers  and  differ- 
ent numbers  of  printers.   This  study  is  a  model  to  cal- 
culate these  costs  and  to  make  a  recommendation  on  an 
optimal  type  and  number  of  printers.   The  costs  to  be 
compared  are  those  incurred  by  the  system  to  include  only 
the  printing  operation  and  the  cost  of  the  queue  that  the 
GARD  system  builds  up. 


II.   APPROACH  TO  THE  PROBLEM 

In  order  to  determine  the  optimal  type  and  number  of 
printers,  the  costs  associated  with  the  various  setups 
must  first  be  determined.   Since  the  message  traffic  flow 
through  the  system  behaves  like  a  waiting  line,  the  appro- 
priate type  of  model  which  best  fits  the  situation  is  a 
queueing  model. 

A.   THE  QUEUEING  MODEL 

To  measure  the  costs  connected  with  a  queue  the 
average  number  of  customers  in  that  queue  must  be  measured, 
In  this  case,  the  customers  are  messages.   It  seems 
reasonable  that  the  relative  cost  of  a  customer  in  a  queue 
is  proportional  to  the  importance  of  the  customer.   This 
importance  of  the  customer  is  measured  in  communications 
by  assigning  a  precedence  or  priority  to  the  message.   A 
flash  message  is  the  most  important,  then  immediate, 
priority,  and  finally  routine.   The  model  must  somehow 
measure  the  average  number  of  messages  of  each  precedence 
in  the  queue.   In  queueing  theory  this  type  of  queue  is 
called  a  priority  queue.   Since  the  alternatives  under 
investigation  include  up  to  a  maximum  of  four  printers, 
the  theoretical  problem  is  thrown  into  the  category  of  a 
multi— channel ,  or  s  server  queueing  model.   Also,  since 
the  arrival  rates  and  service  rates  for  the  messages  of 


different  precedence  are  not  necessarily  equal,  the 
model  selected  from  the  theoretical  library  of  queueing 
models  must  maintain  provisions . for  these  characteris- 
tics.  To  summarize,  the  model  to  bo  used  must  be  a 
multi— channel  priority  queue  with  provisions  for  unequal 
arrival  rates  and  unequal  service  rates  among  the  differ- 
ent precedences.   Because  this  particular  combination  of 
characteristics  makes  the  derivation  extremely  complica- 
ted, the  solution  to  this  problem  is  not  found  in  the 
literature.   For  example,  Cobham  [Ref.  2]  derived  results 
for  a  single  server  priority  queue  with  Poisson  arrivals 
and  a  general  but  unequal  service  time  distribution  for 
the  different  priorities.   He  also  gave  results  for  a 
multi— server  priority  system  with  Poisson  arrivals  and 
the  service  times  exponentially  distributed  with  the  ser- 
vice rate  equal  for  all  priorities. 

The  author  did  some  work  on  a  simpler  model  that  in- 
cluded two  priorities,  one  channel,  Poisson  arrivals,  and 
exponential  service  times.   These  results,  found  in 
Appendix  A,  are  not  new  but  the  approach  is  different. 
Attempts  to  solve  the  same  problem  with  two  servers  could 
not  be  solved  using  the  same  technique.   It  was  obvious 
that  this  type  of  problem  becomes  extremely  cumbersome 
and  uniquely  complicated  whenever  the  number  of  channels 
exceeds  one. 


Since  none  of  the  theoretical  models  fit  the  problem, 
a  computer  simulation  programmed  in  the  General  Purpose 
Simulation  System/360  (GPSS)  language  was  written  and  is 
included.   The  basic  idea  in  the  simulation  is  to  use 
what  the  language  calls  a  storage  location  to  simulate 
the  printing  operation.   The  user  can  then  merely  set  the 
storage  capacity  to  any  integer  to  correspond  to  the  num- 
ber of  printers  being  simulated.   A  transaction  (message) 
trying  to  enter  the  storage  location  is  permitted  to  do 
so  only  if  the  storage  location  is  not  full  (not  all 
printers  are  printing) .   Each  transaction  occupies  one 
unit  of  storage  (utilizes  one  printer) . 

Messages  are  generated  by  the  program  according  to 
the  distribution  specified  by  the  user.   Each  precedence 
is  generated  separately  but  according  to  a  certain 
theoretical  mix.   A  message  is  then  assigned  a  message 
length  according  to  a  distribution  specified  by  the  user. 
Flash  messages  receive  a  special  treatment  at  this  point. 
A  25  word  flash  message  heading  is  added  to  the  message 
length  of  a  flash  message  to  simulate  an  advance  warning 
that  is  placed  at  the  head  of  all  flash  messages.   Next 
the  message  length  of  messages  of  all  precedences  is  con- 
verted to  a  time  that  the  message  will  require  to  be 
printed.   Then  the  message  joins  a  queue  that  is  waiting 
for  access. to  the  printer  or  printers.   The  above  steps 
are  performed  for  each  precedence  of  message.   It  should 
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be  noted  that  all  the  steps  before  the  message  enters  the 
queue  require  no  clock  time  in  the  simulation  but  the 
list  of  events  is  merely  given  to  show  the  reader  the  se- 
quence of  these  events.   Queue  statistics  for  each  prece- 
dence queue  are  maintained  automatically  by  the  GPSS 
program.   At  this  point  all  messages  from  the  four  queues 
vie  for  the  storage  location  (printer  or  printers).   Mes- 
sages are  then  printed  by  precedence.   If  two  messages  of 
different  precedence  are  competing  for  the  printer  at  the 
same  time,  the  one  with  the  higher  precedence  will  gain 
access  when  a  printer  becomes  available.   The  queue  disci- 
pline is  not  preemptive,  that  is  to  say,  a  message  does 
not  knock  a  lower  priority  message  from  the  storage  loca- 
tion (printer  or  printers) .   Immediately  upon  completion 
of  printing,  the  message  leaves  the  queue  in  which  it  was 
originally  a  member.   This  means  that  the  queue  statistics 
that  are  maintained  actually  include  the  time  spent  wait- 
ing to  print,  if  any,  and  the  printing  time.   The  desired 
average  number  of  messages  in  the  printing  queue  or  being 
printed  is  output  as  average  queue  contents  for  the  queue 
of  each  precedence. 

B .   BREAKDOWNS 

A  convenient  way  to  simulate  a  breakdown  or  a  failure 
of  a  printer  would  be  to  merely  introduce  into  the  queue— 
ing  model  a  transaction  (in  this  case  a  breakdown  instead 
of  a  message)  with  a  precedence  higher  than  any  other 
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customer.   The  user  could  specify  the  inter-arrival  time 
distribution  and  the  mean  time  between  arrivals  would  be 
the  Mean  Time  Between  Failures  (MTBF)  of  the  printer. 
The  service  time  distribution  would  be  the  repair  time 
distribution  with  its  mean  time  as  the  Mean  Time  to  Repair 
(MTTR)  of  the  printer.   However,  because  the  MTBF  of  the 
printers  investigated  is  so  long  in  comparison  to  other 
portions  of  the  simulation,  this  procedure  was  not  used. 
The  simulation  uses  a  time  step  routine  and  to  get  accu- 
rate queue  statistics,  some  of  the  printers  require  the 
time  step  to  be  tenths  of  seconds.   To  simulate  an  adequate 
number  of  breakdowns  would  be  a  very  inefficient  use  of 
the  computer  facility. 

Because  the  ratio  of  the  MTBF  to  MTTR  is  at  least  400 
to  one  for  all  alternatives,  the  following  approach  was 
taken  for  the  alternatives  involving  only  one  printer. 
The  number  of  messages  in  a  queue  before,  during,  and 
after  a  breakdown  plotted  versus  time  might  look  something 
like  Figure  1.   Then  let: 

t  =  stochastically  determined  time  to  repair  printer. 

t  =  time  to  restore  queue  to  equilibrium  and  set  to 

t  =  ki  where  k  is  some  positive  real  number  to  in- 
sure that  the  jam  caused  by  the  breakdown  has 
cleared. 

N.,  .  =  averaqe  number  of  messages  of  precedence  j  in  the 

1  j 

system  in  equilibrium  with  one  printer  and  no  breakdov 
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Nn  .  =  average  number  of  messages  of  precedence  j  in  the 
system  during  the  time  t.   This  is  not  an  equili- 
brium condition  since  equilibrium  does  not  exist 
during  a  breakdown. 

A   =  breakdown  rate  or  1/MTBF 
B 

Now  for  one  breakdown,  the  "wasted  message  time"  caused  by 

the  breakdown  is  Nn.t  —  N,  .t.   Because  the  A^  is  so  small 

Oj      lj  B 

the  following  expression  for  the  average  number  of  messages 
of  precedence  j  in  the  system  for  one  printer  including  an 
average  breakdown  is: 

Nxj  =  Nlj  +  AB(N0jt  -  Nljt) 

A  breakdown  package  in  the  program  when  inserted  simulates 
breakdowns  in  succession.   As  determined  by  the  user,  the 
coefficient  k  determines  the  time  to  insure  that  the  jam 
of  messages  caused  by  the  breakdown  has  cleared.   At  this 
time  another  breakdown  is  simulated.   Then  if  n  breakdowns 
are  simulated  in  succession  and  T,  obtained  from  the  com- 
puter output,  is  the  total  time  of  the  simulation  run,  set 
t  =  T/n  and  the  equation  above  becomes 

Sli  =   Nlj  +  VN0jT/n  -  NljT/n)  ■  Nlj  +   ABT/n(NoD   "  Nlj> 

When  i  (i>l)  printers  are  being  simulated,  such  a  dis- 
aster does  not  occur  when  one  printer  fails.   The  queues 
will  merely  tend  to  increase  during  the  breakdown  to  a 
higher  equilibrium  level.   It  should  be  pointed  out  that 
for  those  arrival  rates  of  messages  considered  in  this 
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study,  all  three  models  of  printers  were  capable  of  han- 
dling the  load  with  only  one  printer.   This  is  to  say 
that  an  equilibrium  exists  for  all  models  even  with  just 
one  printer  when  there  are  no  breakdowns. 
Now  if: 
N. .  =  average  number  of  messages  of  precedence  j  in  the 
system  with  i  printers  working  in  equilibrium  and 
with  no  breakdowns. 
AD  =  repair  rate  or  1/MTTR 

Then  A  /(A   +  A  )  is  the  expected  fraction  of  time  that  a 
a       a  k 

printer  will  be  broken  down.  Let  p  =  Ar./(A_  +  A_)  and  if 
Q  .  is  the  average  number  of  messages  in  the  system  when- 
ever there  are  one  or  less  printers  operating  and  P,  = 

i    i— k     k 
(,  )  p    (1— p)   is  the  long  run  probability  that  there  are 

k  printers  out  of  i  working,  then 
i  i 

n.  .  z  j]      p.  n.  .  +  (i  -    y   P,  )  Q,  . 

k=2  k=2 


and  N. .  is  the  average  number  of  messages  of  priority  j 
in  the  system  when  there  are  i  (i>l)  printers  being  simu- 
lated.  Although  Q, .  is  not  readily  computable  it  is 
evident  that: 

Q-,  •  <  N,  . 
ID  ~   ID 

This  is  true  because  Nn  .  was  calculated  while  simulating 

ID 

a  breakdown  of  only  one  printer.   This  implies  that  the 
repair  rate  during  the  breakdown  is  A   or  the  expected 
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repair  time  is  1/A  .  When  looking  at  the  portion  of  time 
that  all  i  printers  are  failed  in  the  calculation  of  Q, . , 
the  repair  rate  is  iX      or  the  expected  repair  time  is  only 

R 

l/iA„.   This  means  that  Q, .  would  be  less  than  N, .  because 
'   R  Ij  lj 

a  breakdown  would  be  repaired  more  quickly  on  the  average. 
These  heuristic  results  complement  the  computer  simu- 
lation to  save  computer  time  that  would  have  been  mostly 
nonproductive  waiting  for  another  breakdown  to  occur. 

C.   COST  MODEL 

One  of  the  most  difficult  tasks  in  any  model  involving 
costs  is  determining  which  costs  are  relevant  and  what 
values  to  put  on  the  costs  that  are  considered  relevant. 
In  this  model  surely  the  costs  involved  with  a  delayed  mes- 
sage are  to  be  included.   The  initial  cost  of  the  printers 
is  a  relevant  cost.   Maintenance  costs  are  another  item 
that  should  be  investigated.   However,  when  the  Navy  pur- 
chases an  item  such  as  a  printer  in  a  large  quantity,  a 
maintenance  agreement  is  usually  included  in  the  purchase 
price  to  include  repair  parts.   Since  maintenance  personnel 
available  for  repair  work  would  probably  be  the  same  for 
all  alternatives,  these  costs  were  omitted.   For  the  pur- 
poses of  this  study,  then,  maintenance  costs  were  not  used 
in  the  computation  of  the  total  cost. 

If  C.  is  the  cost  in  dollars  per  hour  of  delaying  a 
message  of ' precedence  j  and  N. .  is  the  average  number  of 
messages  of  precedence  j  in  the  system  as  derived  in 
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Section  II. B.,  then  total  cost,  TC,  of  a  selected  alterna- 
tive using  i  printers  is  given  by: 

4 
TC.  =  Y]     N. .C .  +  iP/DL 

j=l 

where  P  is  the  price  in  dollars  of  the  printer  involved 

and  DL  is  the  design  life  of  the  printer  in  hours.   If  the 
concept  of  present  value  and  discounting  of  costs  are  used 

the  formula  becomes: 

4 
PV 


r.  =  V  N.  .C.  +  iP 

1   i=l   1D  J  1-e- 


(m)  (DL) 
J 

where  m  is  the  annual  discount  rate  and  set  at  0.10  by  the 
DOD.   The  latter  form  of  the  cost  equation  considering  pres- 
ent value  was  used  in  the  model. 
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III.   SOURCE  OF  DATA  AND  ANALYSIS 

There  were  basically  three  inputs  to  the  model  that 
were  required:   (1)  message  interarrival  time  distribu- 
tion and  rate  of  arrivals,  (2)  the  message  length  distri- 
bution, and  (3)  cost  or  value  of  reducing  the  writer— to— 
reader  time  of  a  Naval  message. 

A.   MESSAGE  INTERARRIVAL  TIME  DISTRIBUTION  AND  RATE 

It  would  be  natural  to  suspect  that  messages  arriv- 
ing at  a  NAVCOMMSTA  do  so  according  to  a  Poisson  process. 
Oelmann  [Ref.  4]  reported  that  this  was  in  fact  true  for 
data  taken  during  the  period  5-8  March  1971  at  NAVCOMMSTA 
San  Francisco.   He  used  the  Chi— square  goodness— of— fit 
test  to  verify  this  result.   Then  it  is  assumed  that  the 
fraction  of  this  traffic  that  is  transmitted  to  a  partic- 
ular ship  is  selected  randomly  from  all  incoming  traffic 
to  the  NAVCOMMSTA.   This  is  to  say,  each  message  that 
comes  into  the  NAVCOMMSTA  has  a  common  probability  of  be- 
ing destined  for  a  particular  ship.   This  assumption 
guarantees  that  the  traffic  transmitted  to  a  particular 
ship  is  still  a  Poisson  process.   The  probability  that 
any  message  is  destined  for  a  particular  ship  must  be  in- 
dependent of  the  destination  of  any  other  message  sent 
before  or  after.   Since  there  was  no  evidence  to  the  con- 
trary, this  assumption  that  messages  destined  for  a 


particular  ship  arrive  at  the  transmitter  according  to  a 
Poisson  process  was  accepted. 

Given  that  the  messages  for  a  certain  ship  arrive  at 
the  transmitter  according  to  a  Poisson  process  does  not 
guarantee  that  the  output  of  the  transmitter  is  a  Poisson 
process,  however.   This  would  be  the  case  if  service  times 
are  distributed  according  to  the  exponential  distribution 
(see  Burke  [Ref.  1]).   In  this  case  the  service  time  (time 
to  transmit)  is  not  exponential  but  is  distributed  accord- 
ing to  the  message  length  (see  Section  III.B.).   It  should 
be  clear,  however,  that  if  the  traffic  is  light,  the  out- 
put should  approach  a  Poisson  process  regardless  of  the 
transmission  time  distribution.   The  fact  that  traffic 
from  four  channels  operating  more  or  less  independently  is 
lumped  at  the  GARD  unit  also  assists  in  making  the  inter- 
arrival  times  of  messages  at  the  GARD  unit  look  like  an 
exponential  distribution. 

A  simplified  transmission  system  consisting  of  four 
independent  channels  was  simulated  as  described  above  and 
the  interarrival  times  at  the  receiver  were  tabulated  and 
a  Chi— square  goodness— of— fit  test  was  accomplished.   These 
results  did  show  that  the  distribution  of  interarrival 
times  at  the  receiver  and  the  GARD  unit  to  be  from  an  ex- 
ponential distribution. 

The  rate  of  message  arrivals  was  determined  from  data 
taken  aboard  the  USS  Jouett,  a  DDG,  in  October,  1969  [Ref. 
5] .   Although  traffic  intensity  varied  somewhat  during  the 
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day,  an  average  rate  of  arrival  was  found  to  be  1/4  mes- 
sage per  minute.   Another  source  of  data  was  that  taken 
by  Oelmann  [Ref.  4]  at  NAVCOMMSTA  San  Francisco  in  March 
1971.   The  rate  indicated  there  was  approximately  1/2 
message  per  minute  to  all  classes  of  ships.   Since  this 
rate  also  varied  considerably  from  day  to  day  and  during 
the  day,  it  v/as  felt  that  a  precise  rate  would  be  hard  to 
determine.   Therefore  the  rate  of  1/4  message  per  minute 
was  selected  and  thought  of  as  a  fairly  heavy  rate  for 
this  class  of  ship.   It  should  be  noted  that  although 
this  may  be  a  fairly  heavy  rate  for  an  average  value,  it 
is  probably  quite  low  for  some  heavy  traffic  periods  that 
all  ships  encounter. 

B.   MESSAGE  LENGTH  DISTRIBUTION 

The  message  length  distribution  was  required  because 
the  length  of  a  message  is  what  determines  that  time  that 
a  message  will  take  to  be  transmitted  or  printed.   As 
stated  earlier,  in  the  queueing  model  this  turns  out  to  be 
the  service  time  of  the  transmitter  or  printer.   It  is 
clear  that  the  service  time  is  scaled  by  the  speed  of  the 
printer  being  utilized  but  the  message  length  alone  deter- 
mines the  service  time  of  any  model  or  type  of  printer. 
The  speed  of  the  transmitter  is  held  constant  at  100  words 
per  minute,  the  speed  at  which  the  Fleet  Broadcast  is 
actually  transmitted. 

The  message  length  distribution  was  found  by  Oelmann 
[Ref.  4]  to  be  a  right  shifted  hyperexponential 
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distribution.   The  source  of  the  data  was  a  sample  of  mes- 
sages taken  at  NAVCOMMSTA  San  Francisco  on  5  October  1971. 
Another  source  of  data  was  the  Naval  Electronics  Labora- 
tory Center  in  its  World  Wide  Traffic  Demand  Study  data 
bank.   Those  data  were  taken  on  8  October  1970  in  various 
parts  of  the  world.   The  determination  of  the  distribution 
as  a  right  shifted  hyperexponential  was  confirmed  by  the 
Chi— square  and  Kolmogorov— Smirnov  (K— S)  goodness— of— fit 
tests.   The  reason  that  the  distribution  is  right  shifted 
can  be  explained  by  the  message  heading  that  appears  on 
all  messages  regardless  of  precedence.   This  heading  con- 
tains a  data— time  group,  an  originator,  and  the  addressee. 
Figure  2  shows  the  cumulative  distribution  function  (CDF) 
for  the  four  precedences  of  messages. 

Both  sets  of  data  were  accepted  by  the  goodness— of— 
fit  as  originating  from  a  right  shifted  hyperexponential 
but  the  data  taken  at  NAVCOMMSTA  San  Francisco  was  used 
to  determine  the  parameters  for  the  distribution  in  the 
simulation. 

C.   COST  INPUTS 

Another  input  to  the  basic  cost  model  was  the  price 
of  time  that  a  message  was  delayed.   In  the  model  it  was 
the  price  of  time  that  a  message  spent  in  the  buffer  be- 
fore being  printed  and  the  time  being  printed.   Washburn 
[Ref.  6]  performed  two  experiments  to  determine  the  price 
of  time  of  a  Navy  message  while  working  at  the  Naval 
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Communications  Command  in  December  1972.   The  most  rele- 
vant experiment  involved  using  the  Delphi  technique  to 
survey  a  group  of  senior  Navy  Communications  personnel. 
The  technique  involves  feedback  to  the  person  completing 
the  questionnaire  in  between  repeated  personal  estimates. 
After  the  first  survey,  all  participants  are  shown  the 
results  of  that  survey  before  being  surveyed  again.   Thus, 
a  member's  estimate  will  be  based  upon  his  knowledge  of 
all  other  participant's  first  response.   The  quantity  of 
interest  in  the  experiment  was  the  price  of  time  for 
Naval  messages.   This  can  be  defined  as  the  amount  of 
money  the  Navy  should  be  willing  to  pay  to  shorten  the 
writer—to— reader  time  of  a  message  by  one  hour.   All  esti- 
mates were  instructed  to  be  of  average  value.   The  results 
of  the  third  and  final  survey  taken  of  ten  senior  communi- 
cations personnel  gave  median  values  of  $1  per  hour  for  a 
routine  message,  $10  per  hour  for  a  priority  message,  $100 
per  hour  for  an  immediate  message,  and  $1000  per  hour  for 
a  flash  message. 

It  should  be  noted  that  although  the  estimates  varied 
widely,  the  median  values  of  price  had  reached  the  final 
values  after  only  two  surveys.   The  third  survey  merely 
lowered  the  variance  within  the  sample. 

D.   MISCELLANEOUS  INPUTS 

Mean  Time  Between  Failures  (MTBF)  and  Mean  Time  To  Re- 
pair (MTTR)  estimates  for  the  three  models  of  printers 
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were  either  obtained  from  manufacturer's  data  or  from  ex- 
perienced technicians  at  the  Naval  Electronics  Laboratory 
Center  at  San  Diego.   Table  I  summarizes  the  values  ob- 
tained and  used  in  the  model.   The  costs  of  the  printers 
were  also  obtained  from  NELC  and  are  included  in  Table  I. 
Department  of  the  Navy  has  set  the  design  life  of  all 
minor  communications  equipment  such  as  printers  at  eight 
years . 

The  fraction  of  messages  by  precedence  that  were  used 
in  the  model  were  1%  flash  messages,  15%  immediate  mes- 
sages, 54%  priority  messages,  and  30%  routine  messages. 
These  figures  were  taken  from  actueil  data  taken  aboard  the 
USS  Jouett  in  October  1969  and  substantiated  by  other 
studies  to  be  meaningful  figures  (see  [Ref.  5]). 

Characteristic 

Model     MTBF      MTTR         Purchase       -  Design 

Price  Life 

Model  28    800  hrs   2  hrs        $  1500  8  yrs 

Diablo     1500  hrs   3  hrs        $  1900  8  yrs 

TT-624     3000  hrs   3  hrs        $15000  8  yrs 

TABLE  I. 
Miscellaneous  Inputs 
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IV.   THE  MODEL 

The  GPSS  program  as  outlined  in  Section  II  was  uti- 
lized to  provide  the  average  queue  contents  statistics 
that  were  required  for  the  cost  model."  When  the  break- 
down package  was  used,  field  A  of  the  TERMINATE  block 
for  all  the  messages  was  set  to  zero  so  that  a  specific 
number  of  breakdowns  would  terminate  the  simulation  run 
rather  than  a  combination  of  breakdowns  and  messages. 
Field  A  of  the  TERMINATE  block  is  the  amount  by  which 
field  A  of  the  START  block  is  decremented  for  each  trans- 
action that  comes  into  the  TERMINATE  block.   The  simula- 
tion terminates  when  field  A  of  the  START  block  is  less 
than  or  equal  to  zero. 

A.   VERIFICATION  OF  THE  MODEL 

Cobham  [Ref.  2]  gives  the  following  results  for  an 

M/G/l  queue  with  r  priorities  of  customers.   Let  S.  be 

the  service  time  of  a  customer  of  precedence  j .   Let 

F.(t)  be  the  CDF  of  the  service  time  distribution  for 
D 

customers  of  the  jth  precedence  and  A.  is  the  arrival 

rate  for  the  customers  of  the  jth  precedence.   Let  r  be 

the  highest  priority  and  1  be  the  lowest  precedence. 

Then 

r  r 

^  =   E  Xa    and  F(t)  =   E   (A./A)F  (t). 
j=l  J  j=l    J     J 
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Then 


dF(t)/dt 


£    (A  /A)f  (t) 
j=l     J     J 


where  f.(t)  is  the  probability  density  function  of  S. 

r 

t2dF(t)  =  £    (A  VA)t2f  (t)dt 
j=l     :       ^ 


/  t2dF(t)  =  /  £   (A./A)t2f  (t)dt 

o  o  j=l    J       J 

r 

00   ? 
=  (1/A)   £   A    /   tTf.(t) 

1=1    J   o      J 


dt 


but    /  t2f.(t)  =  Var(S.)  +  (E(S.))2 
0      D  D         3 

Then  W.  is  the  average  waiting  time  in  the  queue  of  a  mes- 
sage of  precedence  j  and  is  given  by  the  expression: 


W.    = 

3 


(1/2) A    /      t    dF(t) 
o 


j     x 


i-E 


r+l-k 


k=l 


r+l-k 


J  r+l-k 

1-1, 


k=l 


r+l-k 


where  1/y .  is  the  expected  value  of  the  service  time  or 

E(S.)  of  a  customer  of  precedence  j.   Then  the  expected 

time  in  the  system  of  a  customer  with  precedence  j  is 

given  by 

WS.  =  W.  +  E(S  .  ) 
3  3  3 

A  comparison  using  these  results  and  the  results  from 
the  simulation  for  the  case  of  one  Model  28  printer  was 
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performed.   The  results  of  the  comparison  showed  that  the 
expected  waiting  time  in  the  system  from  the  simulation 
was  within  5%  of  the  theoretical  results  for  the  three 
highest  precedences  and  only  slightly  more  for  the  routine 
messages. 

B.  RESULTS 

All  results  from  the  simulation  runs  are  shown  in 
Table  II.   A  95%  confidence  that  all  values  except  the 
breakdown  statistics  are  within  5%  of  their  true  value 
was  attained.   For  the  N~  .  statistics  a  90%  confidence 
that  the  results  were  within  10%  of  their  true  value  was 
attained.   Table  III  shows  the  average  number  of  messages 
in  the  system  to  include  the  provisions  for  breakdown  of 
the  printers.   The  table  also  shows  final  cost  figures 
for  message  delay,  cost  of  printers  and  the  total  cost 
for  each  alternative. 

C.  DISCUSSION  OF  THE  RESULTS 

It  is  evident  from  the  results  that  the  optimum  num- 
ber of  printers  is  two  and  that  the  TT— 624  is  the  most 
economical  model.   Even  if  the  costs  of  a  delayed  message 
are  cut  by  a  factor  of  ten,  the  same  alternative  is  only 
slightly  less  economical  than  two  Diablo  printers.   The 
primary  driving  parameter  in  the  model  is  the  speed  of 
the  printers.   The  slower  printers  waste  too  much  time 
merely  printing  even  when  there  are  adequate  numbers  of 
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them  to  virtually  eliminate  a  queue  in  the  buffer.   At 
times  during  the  day  when  traffic  is  heavier,  the  choice 
of  the  TT— 624  printer  will  be  even  a  better  one  than  dur- 
ing lighter  traffic.   Some  runs  were  made  at  a  rate  of  one 
message  per  minute  and  these  results  substantiate  the 
above  statement.   It  should  be  noted,  however,  that  the 
alternatives  of  one  or  two  Model  28  printers  could  not  be 
simulated  at  this  high  rate  because  it  is  theoretically 
impossible  for  this  model  to  handle  a  one  message  per 
minute  arrival  rate  with  only  one  or  two  servers. 
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V.   CONCLUSIONS 

The  conclusion  that  can  be  drawn  from  this  study  is 
that  two  TT— 624  printers  is  the  most  economic  choice 
from  among  those  alternatives  investigated .   It  should 
be  noted  at  this  point  that  a  lot  of  the  stated  assump- 
tions in  the  model  are  unimportant  or  irrelevant  in  the 
final  analysis.   The  model  points  out  a  rather  common 
sense  result  —  one  high  speed  printer  is  sufficient  but 
that  a  spare  is  required  when  the  primary  printer  breaks 
down. 

An  alternative  that  was  not  investigated  because  of 
programming  difficulties  is  a  mixture  of  two  different 
models  of  printers.   An  example  is  a  high  speed  printer 
to  act  as  the  primary  printer  and  handle  traffic  during 
the  majority  of  the  time  and  a  Model  28,  v/hich  is  inex- 
pensive, to  maintain  a  minimum  queue  during  the  period 
when  the  high  speed  printer  is  broken  down.   This  would 
involve  a  slightly  more  complex  GARD  unit  and  would  in- 
volve extra  repair  parts  so  it  is  not  immediately 
apparent  if  the  advantages  would  outweigh  the  disadvan- 
tages . 

Space  on  a  ship  such  as  a  destroyer  is  also  an  im- 
portant factor  and  has  a  cost  associated  with  it.  The 
most  economic  alternative  involves  the  saving  of  some 
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space  as  compared  to  the  four  Model  28  teletypes  that 
are  presently  installed  in  a  destroyer.   It  was  not  con- 
sidered in  the  study  but  it  is  pointed  out  here  to  show 
another  advantage  of  the  alternative  chosen. 

It  might  be  argued  that  the  Model  28  teletypes  that 
are  currently  in  service  are  free  but  even  if  this  is 
the  case,  the  Model  28  teletype  is  still  inferior  because 
of  its  slow  speed.   It  should  be  pointed  out,  however, 
that  costs  such  as  retraining  costs  for  a  high  speed 
printer  maintenance  course  are  not  included  in  the  study. 
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APPENDIX  A 
SOME  THEORETICAL  RESULTS 

In  the  initial  stages  of  this  study,  the  author  at- 
tempted to  derive  the  theoretical  results  which  were  not 
found  in  the  literature.   A  simple  problem  for  which  the 
results  were  known  was  first  used  to  test  a  solution 
technique.   Given  an  M/M/l  queue  with  two  priorities  of 
customers,  customers  of  the  higher  priority  have  inter- 
arrival  times  distributed  EXP(A„)  and  service  times  dis- 
tributed EXP(y  ).   Lower  priority  customers  have 
interarrival  times  distributed  EXP  (A-,)  and  service  times 
distributed  EXP(y,).   If  the  queue  discipline  is  first— in- 
first— out  within  each  priority  class,  the  higher  priority 
customers  preempt  the  ones  of  lower  priority,  and  P(i,j) 
is  the  steady  state  probability  that  there  are  i  high 
priority  and  j  low  priority  customers  in  the  system,  then 
the  balance  equations  become: 

(X1+A2)P(0,0)    =y2P(l,0)   +y1P(0,l) 

(A1+A2+y1)P(0,  j)=A1P(0#  j-l)+VJ1P(0,j  +  l)+n2P(l,  j)  for  j>l 

(A1+A2+y2)P(i, j)=A1P(i,j-l)+A2P(i-l,j)+y2P(i+l, j ) f or  i>l, 

oo  oo 

j:>0  and  P(i,-1)=0  and   ^       ^      P(i,:j)=1 

i  =  0      j=0 

oo  -1  oo 

Then  let  Q(i,z„)=Xl  P(i/j)z^  and  R(z  ,z  )=Z^  Q(i,z9)z* 
A      j=0         A  X      i=0 

where  0<_z,  ,  z?<L 
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Then  the  balance  equations  can  be  worked  into  the  form 

ylZl(z2~1)P(0,0)  +  [y2Z2(zl~1)~VJlZl(z2""1)  ]Q(0'Z2} 
-A!zlz2+Z2  (X22l+X1z1+y2z1-X2z2-ll2) 

At  this  point  it  would  be  well  to  point  out  that  although 
this  expression  is  still  in  terms  of  the  two  unknowns 
P(0,0)  and  Q(0,z  ),  the  fact  that 

E        D     p(i,j)=i 

i=0    j=0 

(or  equivalently  that  R(l,l)  =  l)has  not  been  used.   After 
this  transformation  to  the  (z,,z„)  space  it  is  now  hoped 
that  by  approaching  the  point  (1,1)  from  different  direc- 
tions that  more  information  about  the  above  equation  can 
be  found  to  solve  for  these  unknowns. 

Then  if    z2  =  f(z,)  -  1+a ( z,-l) +b (z  -1) 2  and  f  (1)  =  1 
f '  (z^  -  a+2b(z  -1)  f'd)  =  a 

f " (z,)  -  2b  f " (1)  =  2b 

Now  let  R(z,z)  =  R(z  ,f(z,))  and  lim  R(z,,f(z  ))  =  -^ 

Zl"1 
so  using  L' Hospital's  rule  yields 

[y  P(0,0)-y  Q(0,1) ]a+u  Q(0,1) 
■■R(l,f(l))  =  — - 


[-A1]a-fy2-\2 


where  Q(0,1)  -  Q(0,f  (z.^  ) 


z1=l 


but  R(l/f(l))=l  for  all  "a".   Therefore  it  is  deduced  that 
Q(0,1)=  l-A2/y   =  l-p2  and  then  P(0,0)  =  l-A-^y  -A2/y2  = 
1— p,—  p„.   Then  if  a  =  (y?— A~)/A-.,  L *  Hospital ' s  rule  can  be 
used  again.   Again  equating  terms  in  the  numerator  and 
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denominator  that  do  not  have  b  as  a  coefficient  because 
R(l/f  (!))--!  for  all  b,  an  expression  is  obtained  for 


Q'  (0,f  (Z]L)) 


z1=l 


It  is  Q"  (0,f  (1))  = 


Now 


RU-^1) 


P1[(1-P2)2+(A2A1)/U2] 
(1-P2) (l-p1-p2) 


y2(z1-l)Q(0/l)      Q{Qfl) 


A2z1+y2z1-A22l-y2   l-p^ 


Then  the  expected  number  of  high  priority  customers  in 


the  system  is  d/dz-,R(z  ,1) 


A 


Zl=l    y2-A2 


which  is  a  well  known  result  from  queueing  theory  consider- 
ing that  the  high  priority  customers  act  as  if  they  are 
alone  in  the  system  and  are  not  influenced  by  the  low 
priority  customers 
Now 


R(l,z2) 


y2P(0,0)-u2Q(0,z2) 


~A2Z2 


Then  the  expected  number  of  low  priority  customers  in  the 
system  is 


af2R(1'z2>' 


^i 

z2=l    A1 


(Q'  (0,l)+P(0/0)-Q(0,l)) 


yl  Ul 

=  ^(Q-  (0,11-pJ  =  ^'(0,1)  -1 

Al  x       Al 

and  inserting  the  expression  for  Q'  (0,1)  that  was  derived 
above  gives  the  result 
d 


dz2R(1'Z2} 


A1p2/y2+p1-p1p2 


Z2  =  1 


(l-p2) (l-p1-p2) 


This  agrees  with  results  found  by  Jaiswal  [Ref.  3]. 
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Using  the  same  technique  as  in  the  above  problem  to 

solve  an  M/M/2  queue  with  the  same  two  classes  of  customers 

as  in  the  above  problem  led  to  a  more  complex  expression 

for  R(z  ,z  )  that  involves  Q(0,z2),  P(0,0),  and  P(0,1). 

(A)Q(0,z2)     +     (B)P(0,0)     +     (C)P(0,1) 
1'    2  X, +A2+2p    — 2y„/z   — A, z   — A„z 

where 

V-.Z,         y,z  2y 

A=(Zl"2+^I  "  ~k      Ul+A2  +  2^1-Al22-  ij> 

+  (X1+X2+2y2  ~  4f"  A1Z2) 

y,z     y  z  2y 

B=(z  -2+— — i-^)  (l-z0)  (A,+A„+  -) 

1    y2z2     y2       2    12    z^ 

2y1z1 
C=(z2-1)  {—r-   +  y2zx  -  2y2  -  y^) 

2 
Again  setting  z„  =  f(z,)  =  l+a(z  — 1)  +  b(z,— 1)   and  using 

L'Hospital's  rule  gives  a  relation  between  P(0,0)  and 

P(1,0)  as  well  as  the  value  of  Q(0,1).   Again  selecting 

"a"  properly  allows  a  second  application  of  L'Hospital's 

rule  but  this  time  no  additional  information  is  found  about 

the  relation  between  P(0,0)  and  P(1,0)  when  the  terms  with 

a  coefficient  of  b  in  the  numerator  and  the  denominator 

are  equated.   Again  Q'  (0,1)  appears  but  its  solution  is  in 

terms  of  P(0,0).   Then  to  take  3—  R(l/z9)  yields  an  expres— 

dz2 
sion  containing  Q'(0,1)  and  P(0,0).   Therefore  the  final 

expression  for  the  expected  number  of  low  priority  custom- 
ers in  the  system  is  still  in  terms  of  the  unknown  quantity 
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P(0,0)  and  cannot  be  solved. 

E (Number  of  low  priority  customers  in  the  system) 


=  —  R  ( 1 ,  z  _ 
dz2      2 


z2=l 


(A)P(0,0)  +  B 
Pi 


where 


A  =  (2+p2) (2p2+2y1+A1+A2)/(]-y1a/y2-2y1a-2y2) 

B  =  2y1/y2-p2-A1/y2-2+2A1-4y2/y1-2ap1+ [2-p2] [p2y1 (2~a) 

+  y1a(y1/y2~2-A1/y2)+A1]  -(2y1)  (2— p-j-p^  (2  +  p2  )  /  (1-y  ]_a/y2 

-2y1a+2y2) 
where  a  =  (2y„— A„)/A, 

One  possible  use  of  an  equation  such  as  this  is  that 
it  is  possible  that  the  quantity  P(0,0)  may  be  easier  to 
measure  in  a  particular  case  than  the  number  of  low  prior- 
ity customers  in  the  system.   Here  this  equation  could  be 
used  to  find  the  expected  value  that  was  desired. 
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